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Abstract 

We propose a new model of leptons and quarks based on the discrete flavor symmetry T', the 
double covering of j44, in which the hierarchies of charged fermion masses and the mildness of 
neutrino masses are responsible for Higgs scalars. After spontaneous breaking of flavor symmetry, 
with the constraint of renormalizability in the Lagrangian, the leptons have mg = and the quarks 
have the Cabibbo-Kobayashi-Maskawa (CKM) mixing angles = 13°, ^23 ~ 0° and = 0°. 
Thus, certain effective dimension-5 operators are introduced, which induce mg ^ and lead the 
quark mixing matrix to the CKM one in form. On the other hand, the neutrino Lagrangian still 
keeps renormalizability. For completeness, we show numerical analysis: in the lepton sector, only 
normal mass hierarchy is permitted within 3(7 experimental bounds with the prediction of both 
large deviations from maximality in the atmospheric mixing angle ^23 and the measured values of 
reactor angle. So, future precise measurements of ^23 1 whether ^23 45° or |023 — 45° | — )• 5°, will 
either exclude or favor our model. Together with it, our model makes predictions for the Dirac CP 
phase, which is almost compatible with the global analysis in la experimental bounds. Moreover, 
we show the effective mass \mee\ measurable in neutrinoless double beta decay to be in the range 
0.04 < |mee|[ey] < 0.11, which can be tested in near future neutrino experiments. 
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I. INTRODUCTION 



In the Standard Model (SM) of particle physics with a single Higgs there are enor- 
mously various hierarchies of quark and lepton Yukawa couplings, that is, rat/in^ = Ut/Uu ^ 
10^, mT-/me = yr/ye ~ 10'^ ctc. In addition, there are several fundamental questions in 
Yukawa sector: why the top quark is uniquely so big compared with the other fermions, 
why neutrino masses are so small and so mild compared with the other charged fermions, 
why both the three leptonic mixing angles and one quark mixing angle are large, while the 
two quark mixing angles are so small. In some sense, our understanding of fermion masses 
and mixing angles remains at a very primitive level. One of the approaches often adopted 
in the understanding for a possible solution for the flavor puzzle consists of the introduction 
of family symmetries which constrain the flavor structure of Yukawa couplings and lead 
to predictions for fermion masses and mixings. We propose a simplified way to address 
those questions in a non-Abelian discrete symmetry T' [ll, by both introducing six types of 
Higgs fields and imposing all Yukawa couplings being of order one. Then the group T' can 
be responsible for the present Pontecorvo-Maki-Nakagawa-Sakata (PMNS) and Cabibbo- 
Kobayashi-Maskawa (CKM) mixing matrices. And the hierarchies of fermions can originate 
from the different Higges. The representations of T' are those of plus three independent 
doublets 2, 2', 2". Similar to ^4 2|-[6| ^, in neutrino sector using four in-equivalent repre- 
sentations 1, 1', 1" and 3 one can obtain the tri-bimaximal (TBM) 7| mixing pattern. In 
the presence of the doublet representations one can naturally describe t 
among the charged fermions and the Cabibbo angle in the CKM matrix 

After the relatively large reactor angle measured in Daya Bay and RENO including 



le mass hierarchy 



Double Chooz, T2K and MINOS experiments [10|, the most recent analysis based on global 
fits ll| of the neutrino oscillations enters into a new phase of precise determination of mixing 
angles and mass squared differences, indicating that the TBM mixing for three flavors should 
be corrected in the lepton sector: their allowed ranges at la (3cr) from global fits are given 
by 



,0+0.44° (+1.30°) ^ _ _ Qnn°+66° (+60°) a _ QQ Qfi°+o.8i° (+2.53°) 

^ -0.46° (-1.47°) ' ' 

= 40.0°;2-i: © 50.4°+} j: la, (35.8° ~ 54.8° Sa) 



(713 — O-UU _o.46° (-1.47°) ' "CP — JUU _-^3go (_3oo°) ' ^^12 — OO.OU _q .^go (_;^ 27°) ' 



^ The finite group ^4 describes the even permutations of four objects and possesses. 
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ir. .70+0.070 (+0.222) 
^•^'-^-COe? (-0.197) ' ^^^^^^ni 
. ,.7+0.042 (+0.185) 
-0.065 (-0.222) ' 

where Amg^j = m\ — ml, Am\^^ = — m\ for the normal mass hierarchy (NMH), and 
Am^tm = I'^s ~ "^2 1 ^he inverted one (IMH). While there are the large values of the solar 
mixing angle 6*801 = ^12? the atmospheric mixing angle 6'atm = 6*23 and the reactor mixing 
angle 6'i.cac = 6*13 in the lepton sector, in the quark sector the Cabibbo angle and the other 



two small quark mixing angles, e.g., at la level 1^ read: 



e\, = (13.03 ± 0.05)° , el, = (2.37^°:°^)° , el, = (0.20!°;°^)° , = (67.17l|H)° .(2) 

The discrepancy of the mixing angles in Eqs. ([1]) and ([2]) may tell us about some new 
flavor symmetries of quarks and leptons. And it is well known that the mass spectrum of 
the charged fermions exhibits a strong hierarchical pattern [see Eq. (!22|) ]. unlike that of 
neutrinos which shows a mild hierarchy. These facts may provide a clue to the nature of 
quark- lepton physics beyond the SM ^. Therefore, it is very important to find a natural 
model that leads to the observed flavor mixing patterns for quarks and leptons. In the 
present article we shall build such a model to emphasize the leptonic mixing parameters and 
the quarks one. 

In this work we propose a new model based on flavor symmetry T' that can accommodate 
quarks and leptons in the same framework invariant under SU{2)l x U{1)y x T' . We intro- 
duce six types of Higgses in the Yukawa sector to depict the mass hierarchies of fermions, and 
simultaneously the mixing parameters of the leptons and those of the quarks. In addition, 
in order to simplify our model and to remove the unwanted Yukawa terms appearing in the 
Lagrangian, we impose a continuous global symmetry which cannot be gauged. After spon- 
taneous U{l)x breaking, to avoid Goldstone bosons it has to be explicitly broken down to a 
subgroup. We stress that in our model CP invariance is originally explicitly broken in the 
Lagrangian level by the complex Clebsch-Gordan coefficients, even though all parameters in 
Lagrangian being real are imposed. The lepton sector in our model can not only naturally 
explain large deviations from the TBM but also provide a possibility for low-energy CP 
violation in neutrino oscillations and the mildness of neutrino masses. 
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The dada ([T]) and ([2|) at la seem to disfavor the maximal mixing in the atmospheric mixing angle 



indicating that it starts to disfavor the Quark-Lepton Complementarity 13| on + ^^23 = 45°, even it is 
not significant yet. 
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The paper is organized as follows. In the next section, we present the particle content 
together with the flavor symmetry of our model and the mass terms of both neutrino and 
charged fermion sectors after flavor symmetry breaking. In Sec. Ill, we show the neutrino 
masses generated in a type seesaw-I and their mixing angles and CP violation as well as the 
CKM matrix. In Sec. IV, for completeness, we show numerical analysis in the lepton sector. 
Then, we give the conclusion in Sec. IV, and briefly mention about the VEV alignments and 
spontaneous CP violation, as an example, in Appendix. 

II. THE MODEL 

In the absence of flavor symmetries, particle masses and mixings are generally undeter- 
mined in a gauge theory. Here, we present a discrete symmetry model based on a T' flavor 
symmetry for leptons and quarks in order to depict the mass hierarchies of charged fermions, 
the mildness of neutrino masses, and simultaneously the present mixing parameters of the 
neutrino oscillation data and those of the quarks. Moreover, we describe the model to un- 
derstand CP violations in the lepton sector which is imperative, if the baryon asymmetry 
of the Universe (BAU) originated from leptogenesis scenario in the seesaw models [l^ . 



Here we recall that T' is the symmetry group of the double tetrahedron [l|,[l5|. The group 
T' has 24 elements and has two kinds of representations. It contains the representations of 
A^: one triplet 3 and three singlets 1, 1' and 1". When working with these representations 
there is no distinction between the group T' and the group A^. In particular, in these 
representations, the elements of T' coincide in pairs and can be described by the same 
matrices that represent the elements in A4. The other representations are three boublets 
2, 2' and 2" . The representations 1', 1" and 2', 2" are complex conjugated to each other. 
Note that /I4 is not a subgroup of T', since the two-dimensional representations can not 
be decomposed into representations of ^4. The generators 5* and T satisfy the relation 
S"^ = R,T^ = [STY = = 1, where i? = 1 in case of the odd- dimensional representation 
and R = —1 for 2, 2', 2" such that R commutes with all elements of the group. In the 
three-dimensional unitary representation, there are abelian subgroups of T' : Z^, Z4 and Zq 
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symmetries, which are generated by the elements 



T 



/ 1 

a; 




TST^ 



\ 



( 
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-1 2uj 2u;2 
2a;2 _i 2uj 
2uj -1 2uj^ 



\ 



(3) 



-12 2 
2-12 
\2 2 -IJ y 

respectively, where to = e*^'^/^ is a complex cubic-root of unity. Especially, the elements 
T and TST"^ are of importance for the structure of our model. The group T' has seven 
irreducible representations, one triplet 3, three doublets 2, 2', 2" and three singlets 1, 1', 1" 
with the multiphcation rules V ® = ® V = r-'^^' for r = 1, 2, V ^3 = 3 ^ V = 3, 
2^' ® 2'= = 3 © 1^'+'= = rJ■+^ 2^' ® 3 = 3 ® 2^' = 2 © 2' © 2" and 3 ® 3 = 3, © 3^ © 1 © 1' © 1", 
where j, /c = 0, ±1 and we have denoted 1° = 1,1^ = 1', 1^^ = 1" and similarly for the 
doublet representations. The sum j + k is modulo 3. The Clebsch-Gordan coefficients for 
the decomposition of the product representations are shown in Ref. [15 



16|. 



We extend the standard model (SM) by the inclusion of an T'-triplet of right-handed 
S'f/(2)L-singlet Majorana neutrinos A''^, and the introduction of six types of scalar Higgs 
fields: the SM S'f/(2)i-doublet Higgs bosons $, which we take to be T'-triplets 3 repre- 
sentation, the other two S'f/(2)i-doublet Higgs bosons H,G, which are distinguished from 
$, \1/ by being T'-doublets 2 representation, and another 5'[/(2)£,-doublet 2 of Higgs boson 
7], which is a T'-singlet 1 representation, finally an S'f/(2)i-singlet T'-triplet 3 Higgs field x'- 




Hk 



Gk 




Xji 



V 




(4) 



where j = 1,2,3 and = 1,2. According to the above Higgs scalars, we impose T' fiavor 
symmetry for leptons and quarks. And, we levy all Yukawa couplings are of order one which 
implies that all the hierarchies of fermions appearing in the Lagrangian are responsible 
for the Higgses we introduced economically. In addition, after spontaneous breaking of 
fiavor symmetry, the VEVs (vacuum expectation values) of such fields need spontaneous CP 
violation. The representations of the field content of the model under SU{2) x f/(l) x T' 
are summarized in Table-Hl where each fiavor of lepton doublets is assigned to one of the 
three T'-singlet representations: the electron-fiavor to the 1, the muon fiavor to the 1", 
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and the tau flavor to the 1', and Ql denotes left handed quark SU{2)l doublet and tr,£r 
and tR^lAn (fe/j,!^/?) are the respective SM right handed lepton and u-type (rf-type) SU{2)l 
singlets, respectively. Here, the down-type fermions Vr and the up-type fermions Ur 
are assigned as T'-doublets and right handed gauge singlets: 

U„ [ . (5) 

m \^ Sr \^ Cr 

In the presence of three T'-triplet Higgs scalars $, \E', x and two T'-doublets Higgs scalars 
H, G, Higgs potential Lagrangian involving interaction terms among $, x, H, G, \1/ and 
X, which would be written as V{^x), V{^H), V{^G), V{^H), V{^G), V{Hx), V{Gx) and 
y(\E'x), would be problematic for vacuum stability. Such stability problems can be naturally 
solved, for instance, in the presence of extra dimensions or in supersymmetric dynamical 
completions In these cases, those interaction terms are either disallowed or highly 
suppressed. In our model, the T' flavor symmetry is spontaneously broken by those T'- 
triplet and doublet scalars and T'-singlet scalar. From the condition of the global minima 
of the scalar potential, we can obtain vacuum alignments of the flelds x^^y"^ j H,G relevant 
to achieve our goal. The Higgs potential of our model contains many terms, which is listed 
in Appendix El Eqs. (lA3l) -( lA9l) . We spontaneously break the T' flavor symmetry by giving 
non-zero vacuum expectation values to some components of the T'-triplets ^ and $. We 
take the T' symmetry breaking scale to be above the electroweak scale in our model, i.e., 
(x) > {^^)- As seen in Appendix |Al the minimization of our scalar potential gives the 
following vacuum expectation values (VEVs): 



(X) = ^,e^ni, 1, 1) , (<f °) = ^(1, 0, 0) , = ^(1, 0, 0) , 

The SM VEV v = (V2G^)-i/2 ^ 246 GeV results from the combination = + 
^1 + vjj^ + Vq^ + f|) where k = 1,2. The non-zero expectation value (x) ~ (1)1)1) 
breaks T' symmetry down to a residual symmetry which is generated by the group 
element TST"^. The non-zero expectation values {^^) ~ (1,0,0) and ~ (1,0,0) break 
T' symmetry down to its subgroup Z3 which is generated by the group element T. The 
non-zero expectation value (77°) = Vn/V2 does not break the T' symmetry, because it is 
T'-flavorless. The non-zero expectation values {H^) and {G^) break T' — )■ nothing with 
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TABLE I: Representations of the fields under T' and SU{2)l x J7(l)y. 



Field Lr,L^,Le Ql tr,£r IrMr bn^Vn Nr x V ^ H ^ G 
T' 1', 1", 1 3 1', 2" 1', 2' 1', 2" 3 3 1 3 2 3 2 

S[/(2U X [/(l)y (2,-i) (2,i) (1,-1) (2,1) (2,-i) (1,0) (1,0) (2,i) (2, i) (2,i) (2, i) (2,i) 



hierarchical breakings. In addition to T' flavor symmetry, we impose an additional symmetry 
U{l)x which is continuous global symmetry, where Le,Lfj,,Lr, Tr, Sr, Br and Vr carry 
X = 1 and $, G carry X = —1, while all other fields have X = 0. So this non-fiavor 
symmetry forbids some irrelevant the SU(2)l x f/(l)y x T' invariant Yukawa terms from 
the Lagrangian (see later). And Since Goldstone bosons resulting from spontaneous U{l)x 
breaking via ($), (G) ^ are not allowed phenomenologically, so the additional symmetry 
U{l)x has to be explicitly broken ^ down to a subgroup Z2 under which $ — )■ — $, G — )■ —G 
and the fields Lg, L^, Lr, tr, Sr, Br and Vr also switch sign. 

In our Lagrangian, we assume that there is a cutoff scale A, above which there exists 
unknown physics. 

A. The neutrino sector 

The Yukawa interactions {d < 5) in the neutrino sector invariant under SU (2)^ x U (l)y x 
T' can be written as 

- = yr^.($iV^)i' + y'2L,mR)r' + y'M^NR)^ 

+ ^M(N^Nr), + lyUM.NR)sa + h.c, (7) 

where $ = ir2$* and T2 is a Pauli matrix. Note here that there are no dimension-5 opera- 
tors driven by x fi^ld in the neutrino sector, and the above Lagrangian in neutrino sector is 
renormalizable. In this Lagrangian, each flavor of neutrinos has its own independent Yukawa 
term, since they belong to different singlet representations 1', 1", and 1 of T': the neutrino 
Yukawa terms involve the T'-triplets $ and Nr, which combine into the appropriate singlet 
representation. The right-handed neutrinos have an additional Yukawa term that involves 

^ In Appendix El there are interaction terms {G^ H){G^ H) which break explicitly U{l)x to 

remove the unwanted Goldstone bosons in the low energy spectrum. 
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the T'-triplet SM-singlet Higgs x- The mass term |M(A'"^A^r)i for the right-handed neutri- 
nos is necessary to implement the seesaw mechanism by making the right-handed neutrino 
mass parameter M large. The additional symmetry U{1)r, as shown before, guarantees 
that the SU(2)l x f/(l)y x T' invariant Yukawa terms Le^^^r^Nn are forbidden from the 
Lagrangian. 

After the breaking of the flavor and electroweak symmetries, with the VEV alignments 
as in Eq. (jH]) the Dirac neutrino and right-handed neutrino mass terms from the Lagrangian 
d?]) result in 



- C 
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V2 
M 



- Nl,Nj^ - Nf,Nn2 - WiNr2 - W2NR1 - Wi^Rs - Ws^ri } + h.c. . (8) 

Then, the neutrino Dirac mass terms and the right-handed Majorana neutrino mass terms 
are expressed as 



rriD 



v<s,e 



47 



^y- ^ 

y:^ 



V2 



(9) 



ivRVxe''' 



-IvRVxe''' 
IvRVxe''' 



M 



\ 



lip 



■sVR^xe 

where 2/123 and M, y'^, are real positive variables 



ivRVxe''' 
- bn^xe''' 
IvRVxe'^ 



(10) 



B. Charged fermion sector 

In the charged fermion sector, the Yukawa interactions {d < 5) including dimension-5 
operators driven by the x field, invariant under SU{2)l x U{1)y x T', are given by 

■^Yuk ~ -^Yuk + -^Yuk + -^Yuk ; (H) 
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where 

a(s) Yd Y'^(s) 

+ ^[(QL^)3X]i"&iJ + -fQL{HVR)r,x + ^Ql{HVr)^x + h.c, (12) 
= yt{QL^)v'tR + Y^QL{GUR)^ 

a(s) Yu y'^M 

+ ^[{QL^hx]i"tR + -fQL{GUR)rX + -^Ql{GUr)sX + h-c , (13) 
-^Yuk = VrLrri TR + Y^L^{H£R)r> 

+ ^-L^[{H£r)^x]v' + ^Ie[(i/£'R)3X]i + h-c, (14) 

with G = iT2G*. In the charged- lepton sector tau lepton involves the T'-singlet r] and the 
T'-singlet right-handed charged-lepton tr, while muon lepton has the T'-doublet H and the 
T'-doublet right-handed charged-lepton Sr and there is no corresponding electron lepton 
term in renormalizable terms, indicating directly electron mass can be generated by the 
dimension-5 operators driven by x field. And the r-mass is generated upon the breaking of 
T'— invariant. Thus, the third family in charged leptons is different from the two, muon and 
electron. On the other hand, in the quark sectors the use of T'-doublets H, G and triplets 
$, \& Higgses can allow the third family to differ from the first two, and thus make plausible 
the mass hierarchies rrit ^ rric ^ rriu and nif, ^ ^ j^. And in the renormalizable 
terms the 6-quark and t-quark masses are generated upon the breaking of T' — )■ Z3. Mass 
terms of the quarks have two independent Yukawa terms with different couplings {yb, Yd) 
and {yt, Y^) for down-type quark and up-type quark, respectively, all involving T'-triplet 
Higgs fields $, \1/ and doublet fields H, G. The T'-triplet $ is shared by both the three 
neutrino Yukawa terms and top-quark Yukawa term. The T'-doublet H is involved by the 
terms associated to the T'-doublet right-handed lepton Sr and down-type quark Vr. As 
mentioned before, the above Yukawa Lagrangian has the additional symmetry U{l)x- This 
non-fiavor symmetry (continuous Global symmetry) ensures that the SU(2) x t/(l) x T' 
invariant Yukawa terms Ql{GVr)s,Ql{GSr)s, Ql{HUr)^, {QL^)i"bR and {QL^)v'tR are 
absent from the Lagrangian. 

In the charged fermion sectors from the Lagrangian ffTTl) . after the breaking of the fiavor 
and electroweak symmetries, with the VEV alignments as in Eq. ([6]) the up-type quark and 
down-type quark mass terms result in 



where 



= yb-j=bLbR + -j=^iVH,dLdR H —sl [(IrVh^ + srVh^) + vh^OlSr 

V ( — — — — 

+ -^^mj^ bibR + 5^6^ + dibR + c^L^i? + m^^ slSr + ^ls^j 



+ dic/R + 77121 ^LdR + feirfi? j + h.c. 



(16) 



1 - i- 



= yt-j=tLtR + -^jt^Ga^iCij + Wg.ClUr + — (^^R^Ga + CrVg^ 

V ( — — 

H ti^tR + CLtR + m^g Uit/j + UlCr + m22 ClCr + tiCR 



V2A l 

+ ulUr + m2i clUr + tin/?^ [> + h.c. 



V Y _ _ 



(17) 



(18) 



with v^, = f^e*'',ti$ = f$e*'>',-UGfc = vc^e^"^ , vh^ = vu^e^''^ {k = 1,2) and = v^e^'-^ . In the 
above Eqs. ffT^ . f lT7|) and f lTS]) . the entries m'^-,m\ - and mf^- are given in Appendix |Bl Then, 
the down- type quark mass matrix A^^ is given by 

( lYdVH, ^ 

1 

\^ YdVH2 VbV^ j 

l^/Diag(m,,m„m,)V^^^ . (19) 



lYd^H, 
^YaVH, ^-^Y^vh, 



tip 



"^11 "li2 ?Tii3 



^21 ^^22 ''^23 
\ "^fl "^^2 "^fs / 



And, the up-type quark mass matrix hA^ can be exphcitly expressed as 





lY^VG, 
V2'Diag(m„,mc,mt)V^"^ . 



tip 



^ m\i ^ 



V2A 



m*2i m*22 m*23 



\ "^31 "^32 "^*33 / 



(20) 



Finally, with the VEV alignment in Eq. ([6]) the charged- lepton mass matrix Aii can be 
explicitly expressed as 

'^O o\ / mi n\ 

T2 -^^-^^^^^ +7il 





Yy,vu2 Yf,VHi 

y y^t;^ y 
Vl Diag(me, m/^m^) 



v 










m2i 


p 















(21) 
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In Eqs. (US]), ([20D and dnD, V[ and are the diagonalization matrices for Mf. 

There exist several empirical fermion mass hierarchies in the up- and down-type quark 
and charged- lepton sectors calculated from the measured values jsj : 



m„ = 2.4 MeV = 1.27 GeV rrit = 171.2 GeV 

TUd = 4.8 MeV TUs = 104 MeV rrib = 4.2 GeV 
me = 0.511 MeV = 105.7 MeV = 1.777 GeV (22) 



which implies that the possible quark-lepton symmetry [17|] is broken by the masses of quarks 
and leptons. Thus, it is not expected that the known quark mixing pattern is transmitted to 
the lepton sector in the exactly same form. In addition, a key point inferred from Eq. is 
that the mass spectrum of the charged leptons exhibits a similar hierarchical pattern to that 
of the down-type quarks except for the electron mass which is much smaller than d-quark 
one, unlike that of the up-type quarks which shows a much stronger hierarchical pattern 
and top-quark is uniquely biggest. Further, there is another interesting empirical relation 



1 , , 1 

2 



« ('^ V « 3 f- ) . (23) 



which has been known for quite a long time 18|. For instance, in terms of the Cabbibo angle 



A = sin^c ~ \Vus\, the fermion masses are scaled as {meiim^) ^ (A^,A^) rrir, {md^nis) ^ 
(A^, A^) nib and {muiTnc) ~ (A®, A^) m-f, which may represent the followings: (i) there is 
at least one Higgs scalar shared by both charged-lepton and down-type quark sectors, or 
(ii) the mixing matrix of the charged lepton sector is similar to that of the down-type quark 
sector, and (iii) the CKM matrix is mainly generated by the mixing matrix of the down-type 
quark sector. 

One of most interesting features observed by experiments on the charged fermions is that 
the mass spectra of quarks and charged leptons are strongly hierarchical, i.e., the masses 
of third generation fermions are much heavier than those of the first and second generation 
fermions. For the elements oi Aif given in Eqs. (l20l) and (fT9l) . taking into account the most 
natural case that the charged fermion masses have the strong hierarchy nit ^ "^c ^ '^u 
{mb{T) ^ ms(^fj_) ^ nid{e)) as well as Eq. (j22l) . we make a plausible assumption 

Vtv^ > VbV^ = VrVr^ = YuVg^ > YdVn^ > YuVg2 = YdVH^ ■ (24) 

Then and can be determined by diagonalizing the matrices M-jM^^ and M^^Mf, 
respectively, indicated from Eqs. (|2T|) . fl20|) and (fT9|l . Especially, the mixing matrix V/ 
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becomes one of the matrices composing the PMNS and CKM ones and it will be shown 
later. A general diagonalizing matrix V[ can be parameterized in terms of three mixing 
angles and six phases: 



V 



v 



7 > 



(25) 



where s{ = sm9{, c{ = cos9{ and a diagonal phase matrix Pf = diag(e*^i , e*^2 , e*^3 ) which 
can be rotated away by the phase redefinition of left-charged fermion fields. 



1. The down-type quark sector and its mixing matrix 

First, we consider the down-type quark sector. From Eq. fll9p we see that the down-type 
quark mass matrix Aid can be diagonalized in the mass basis by a biunitary transformation, 
V;^^MdVi^'^ = Diag( ■md,ms,mb). The matrices and can be determined by diagonal- 
izing the matrices Aid-Mli and Ail^Aid, respectively. Especially, the left-handed down-type 
quark mixing matrix becomes one of the matrices composing the CKM matrix such as 
VcKM = yl^yt [see Eq. (gl])]. From Eq. ^ the hermitian square of the mass matrix for 
down-type quark JUdMd can be obtained. And, from Eqs. (IT^ and the mixing angles 
and phases can be expressed in terms of Eq. (12^ as 

+ 0^^^ + larg(.-l) , (26) 

where the parameters and are positive real numbers of order unity. The empirical 
relation Eq. fl2^ for down-type quark can be satisfied by setting as follows 

A \2yb 3ybJ A 3yb V2 

where the parameters Ad and Bd are positive real number of order unity. Note that the third 
relation in Eq. f l27|) comes from the renormalizable terms. Then, from the above relation 
the mass squared eigenvalues are written as 

V^^MdM^Vi ^ Diag {ml ml ml) 

^Diag(^^ri^{risin(pi2-^) + 1^2Cos(pi2-(p)}, ^y^l) , (28) 

12 



y a y s 

+ TT^- So, one can obtain the 



where FrfW^a = V2X^ybVq;, and Yi = ^ + ^ + Y^, 
measured value of md/mi, ~ A^. Then, we can obtain the mixing matrix of the down-type 
quarks: under the constraint of unitarity up to (9(A^), it can be written as 



1 - 



A2 



2 



1-^ 

2 



-i(f+f) 



fidA^e^2 
AdA^e^d+f) 
1 



Pd + O(A^) . (29) 



2. The up-type quark sector and its mixing matrix 



Next, let us consider the up-type quark sector to obtain the realistic CKM matrix. From 
Eq. (I2U]) and Eq. (I24p the hermitian square of the mass matrix for up-type quark A^^A^]^, 
with the condition given in Eq. (12^ the mass squared eigenvalues are written in a good 
approximation as 



2 ^\^2 2 



2 -'--1^2 2 



(30) 



Then, one can set mu/rrit ~ vq^I"^^ — rrLc/mt ~ Vd/vq, ^ A^ for equal amounts of 

Yukawa couplings. And, the mixing angles and phases can be expressed in terms of Eq. (!24|l 

as 



' A V2yt ^yt 



2 



vr 

2 ' 



^2 — 



A 3yi 
" 2 ' 



0" ~ ^ 



A V2F, 



3K 



- + ^ 
2 2 



(31) 



Due to v^J K ~ A^ in Eq. ( 127|) and the measured value of mu/mc ~ vq^I'^Gx ~ -^^ in Eq. (!22|) . 
it is impossible to generate the Cabbibo angle, A ~ |Kjs|, from the mixing between the first 
and second generations in the up-type quark sector: if one sets 1(^/^)12! = ^3, then from 
Eqs. (!30!) and (13T1) one obtains |(V2')i2| ~ "^xM ~ -^^5 discrepancy with the measured 
A ^ And similar to Eq. fl27|) . from Eq. fl3T]) one can set 



^3" - Cuy? 



(32) 



where Au,Bu, are positive real numbers of order unity. Then, plugging Eqs. fl3Tl) and 
(!32|) into Eq. (l25ll the up-type quark mixing matrix can be written, under the constraint 



13 



of unitarity up to 0{\^), as 



1 



1 



tyg + TT 



Pu + 0(A' 



(33) 



which indicates that the mixing matrix in the up-type quark sector can affect, at most, the 
next leading order contributions in A. Flavor-changing neutral currents (FCNCs) in up-type 
quarks may restrict Au,Bu,Cu < A, which we will not discuss. 



3. The charged lepton sector and its mixing matrix 

Similar to the quark sector, from Eq. (12T!) we see that the charged lepton mass matrix 
A4i can be diagonalized in the mass basis by a biunitary transformation, V^^A^^V^^^^ = 
Diag(me, m^, rrir). The matrices Vl and can be determined by diagonalizing the matrices 
AiiAil and Ai^Aie, respectively. Especially, the charged lepton mixing matrix becomes 
one of the matrices composing the PMNS matrix such as f/pMNS = ^l^v^ respectively [see 
Eq. (!60!) ]. From Eq. (12T|) the hermitian square of the mass matrix for charged lepton J\A.(,M\ 
can be obtained. From Eqs. fl2T|) and fl25l) the mixing angles and phases can be expressed 
in terms of Eq. (12^ as 

^',^4a.g(..-l), (34) 

and </)f = 02 = 0, 9{ = 62 = 0, where the parameters are positive real numbers of 

order unity. Then, from the above relation the mass squared eigenvalues are written, in a 
good approximation, as 

V^f Al,>llV^(^) ^ Diag K, m% m?) ^ Diag Qf^ (^)', ^F^i^, \yy^ .(35) 

Due to the ratio of the measured values me/m^ ^ A'^ and m^/mj. ~ A^, they can be expressed 
in terms of Eq. f l2T]) as mf.jm^ ~ Y^v^j^Y^bC) and m^/m^ ~ Y^Vh2/ {UrVr))- And, one can 
express 

= A,X' , (36) 

where the parameter is positive real number of order unity. Then, we can obtain the 
mixing matrix ¥[ of the charged leptons: under the constraint of unitarity up to 0{X^), it 
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can be written as 



V, 



v 







1 




Pi + 0{X' 



(37) 



It indicates that the effect of mixing in the charged-lepton sector to the PMNS matrix is at 
least less than A'^ ~ 0.2°, and its contribution to the PMNS matrix is negligible because of 
the relatively large reactor angle ^13 measured in Daya Bay and RENO experiments 9|. 



III. PMNS AND CKM MIXING MATRICES AND MASS SPECTRA 
A. Quark Sector 

In the weak eigenstate basis, the mass terms in Eqs. ( IT6|) and (ITTj) and the charged gauge 
interactions can be written as 

- = qtMuq'-R + ^lMmI + qirqi + h-c . (38) 

Let us first consider the quark sector. From Eq. (l38i) . to diagonalize the up- and down-type 
quark mass matrices such that 

yPMfVl = Diag(m;,,m;„m^3) , (39) 

we can rotate the quark fields from the weak eigenstates to the mass eigenstates: 

■u(d) ^Md-)i u{d) u(d) ^Md)iu(d) / 

VL ~^ HL 1 HR ~^ HR ■ V^^) 

Then, from the charged current terms in Eq. f l38p . we obtain the CKM matrix 

VcKM = Vl^V^ . (41) 

From Eqs. (l29l) and ( l33l) we can obtain directly V^'^V^, and by recasting the result with 
the transformations d de*^', c — )• ce*('^3-5+?')^ g _i. se*(</'3-5)^ i te''^'t>i+%+^-i) 
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b — 5> &e*^''^3+2 + 2 f), we can rewrite the CKM matrix as ^ 



CKM 



1-f + C,A3sin(0f-f) 
-A + CuX^ sm{(f)i - f ) 

O(A^) . 



A - CuX^ sin(0^ - f ) 5A3e^^(<^3+f -0 



l-4^ + aA3 sinfc/.^-^: 



-AA^ 



1 



(42) 



where ^ A = - S = - B„, ^' ~ C^A^ cos(f - 0f) and ^ ^ C^A cos(f 



bf). If one 



y a y s 



set — 7- which can be reahzed by Cu - 2y ' 3Y 



A" with n > 1, and 



(43) 



then one can obtain the reahstic CKM matrix in the Wolfenstein parametrization 19|] given 
by 



CKM 



^ 2 

-A 



A AX^{p + ir]) 



1-f AX' 



yAX^l- p + iT]) -AX' 



+ 0{X' 



(44) 



As reported in Ref. [12[ the best-fit values of the parameters A, A, p, t] with la errors are 



A = sin^c = 0.22543 ± 0.00077 
p = 0.144 ±0.025 , 



A = O.SUt^ol^ , 



(45) 



where p = p(l — A^/2) and t] = 77(1 — A^ / 2) . The effects caused by CP violation are always 



proportional to the Jarlskog invariant 20 1, defined as J^p"^^ = lm[VudVcsV*^V*^] ^ A'X^ 
whose value is 2.96^0:17 x 10"^ at la level [12|. 

From Eqs. (128|) and ( 130|) . the observed quark masses respect the following relation 



mrf : nis : nib ^ VHoYd 



Xv^ 



1. 



1 



{^1 sin(pi2 -(p) + Y2 cos(pi2 - ^)} : ^YdVH2 ■ -^V^^^ 



(46) 



^ Note here that since the matrix in Eq. (|42|) has the phase ip dependence which is from x field, when 
Cu ^ A" with n> \ [n: integer) one can re-parameterize and obtain the CKM matrix. 

^ Taking into account the FCNCs in the up- type quarks we may approximate Au,Bu,Cu ^ [see below 
Eq. (133] . 
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B. Lepton Sector 



The mass terms in Eqs. (|8]) and (fTSi) and the charged gauge interactions in the weak 
eigenstate basis can be written in (block) matrix form as, using N'^mo'^l = T^m'^Nfi, 

(47) 



mW 



\n^^MrNr + UZmoNR + liMdR + W-£l7^z/l + h.c. 




9 



+ ^LMdR + ^W^^hY^L + h.c. (48) 

Here £ = (e, r), z/ = (z/e, z/^, z/^), iVfl = {Nri, Nr2, Nrs). 

To find the neutrino masses and mixing matrix we are to diagonahze the 6x6 matrix 



rriD 
ml Mr 



(49) 



We start by diagonahzing Mr. For this purpose, we perform a basis rotation Nr = Uj^Nr, 
so that the right-handed Majorana mass matrix Mr becomes a diagonal matrix Mr with 
real and positive mass eigenvalues Mi = aM, M2 = M and M3 = bM, 

\ 



Mr = U^^MrUr = MUl 



1 



lip 



\ 



^Ke'^ 1 - i 



lip 



^ aM ^ 



U. 



R 



M 
bM 



(,50) 



where k = y'^v^/M. We find a = + + 2k cos (p, b 
diagonahzing matrix 



1 + k' 



2k cos and a 



f \ / e^^ ^ 



R 



3 

1 1 
]_ J_ J_ 



1 

V2 



1 

e'^ 



(51) 



with phases 



^1 



tan 



-Ksmf/) 
I- K cos ip 



and '?/;2 = tan 



KsmLp 



(52) 



. 1 — K cos (p / 

As the magnitude of k defined in Eq. ( l50l) decreases, the phases ipi^2 go to or vr. With the 
basis rotation A'^^ — )■ Uj^NR, the Dirac neutrino mass matrix gets modified to 



mo = moUR 



V2 



1 
2/2 

V y 



(53) 
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where yi = yl/y'^,y2 = ?/2/z/3- At this point, 



-C 



1 



hMdR + -^W;hl^UL + h.c. 



(54) 



mo 
ml Mr I \Nr^ 

Now we take the hmit of large M (seesaw mechanism) and focus on the mass matrix of the 
hght neutrinos My, 

-Cmw = ^i^Myul + hMeiR + -^W-TI-f^UL + h.c. + terms in Nr (55) 

with 



Mu = —mo M^^ m^- 



(56) 



We perform basis rotations from weak to mass eigenstates in the leptonic sector, 



= Vl^iL , Ir = VSiR , 



(57) 



where U^, ^l{r) are unitary matrices chosen so as the matrices 



My = UlMyU: 



-UlmnURM^\UlmnUR) 



Ml = V'^MiVi 



R 



(58) 



are diagonal. From Eqs. (1531) and ([3SD the observed charged lepton masses respect 



me : m^ : mr — 



(59) 



And from the charged current term in Eq. 

as 



we obtain the lepton mixing matrix f/pMNS 



PMNS — l^L '^i' 



Vf'Uy. 



(60) 



The matrix t/pMNS can be written in terms of three mixing angles and three CP-odd phases 
(one for the Dirac neutrinos and two for the Majorana neutrinos) as 8| 



PMNS 



v 



C13C12 



C13S12 



sise 



-i&CP 



"C23S12 — S23Ci2Sl3e*'^^^ C23C12 — S23Si2Si3e*^'^-f' S23C13 



■523'5l2 ~ C23Ci2Si3e 



«0CP _ 



■523^12 ~ C23Si2Si3e* C23C13 



Pu , (61) 



where Py = Diag(e e *'^2/2^ j^-j^ g,, ^ ^i^Q.. and Qj = cos6'i 
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After seesawing, the light neutrino mass matrix is given by 



(62) 



where we have defined an overall scale tuq = v\y'^'^ / {QM) for the light neutrino masses. 
And the overall phase can be rotated away by redefining the light neutrino field. The mass 
matrix A^,y is diagonalized by the mixing matrix Uy, 

Mu = Diag(mj.,,m,,2,mj,3) [/J. (63) 

Here rrij^. {i = 1,2,3) are the light neutrino masses. Interestingly, the mixing matrix Ur in 
Eq. fl5T]) refiects an exact TBM. Therefore Eq. f l62|) directly indicates that there could be 
deviations from the exact TBM if the Dirac neutrino Yukawa couplings do not have the same 



magnitude. In the limit y'( = the mass matrix in Eq. flB2]) acquires a yU-r symmetry 21| 
that leads to = and 6^23 = — 7r/4. Moreover, in the limit yi=y2= Vs (z/i, Z/2 — ^ 1), the 
mass matrix ( 162|) gives the TBM angles and the corresponding mass eigenvalues: 

= , 91^ = 0, 9^, = sin-i (^-L^ , 

3mo 3mo . . 

mi,^ = , 771^2 = 3"^o , = —r- ■ (64) 

a h 

These mass eigenvalues are disconnected from the mixing angles. The neutrino texture in 
Eq. ( I62l) provides naturally the mildness of neutrino masses, because the components giving 
the TBM are multiplied and constrained by neutrino Yukawa couplings. 

Due to in general 1/1,1/2 7^ 1? there are deviations from their TBM values. Moreover, recent 
neutrino data, i.e. 6'i3 7^ 0, require deviations of yi^2 from unity because the contribution of 
V[ in Eq. ( l37|) from the charged-lepton sector is expected to be small, leading to a possibility 
to search for CP violation in neutrino oscillation experiments. To diagonalize the above mass 
matrix Eq. (162|) . we consider the hermitian matrix A^j^A^J,, from which we obtain the masses 



^ The neutrino mass matrix form given by Eq. (|62l) is different from tlie one given in Ref. Q due to the 
opposite sign in f 
numerical results 



opposite sign in front of ^^^5— in 2-3 sector of the light neutrino mass matrix. So, it makes a difference in 
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and mixing angles: 



A B C 
B* F |G|e^'^i 
yc* \G\e-'^'i K j 



Uy Diag(m^^,m^2,m^J Ul 



(65) 



where 



A = p + q + 2gi , F = yl(^p+'^ + r - gi- , K = yl [p + + r - gi + g^ 



B = yAv--^ + ^ 



G 

with 



2/12/2 [V 



4 



93 ■ 
r- gi 



.3 (5-4 + 5-5) 



2a 



C 



yi P 



1,9i 
2 2 



93 + i 



3(^4 - ^5) 
2a 



(66) 



^1 
^3 
9b 



rrin 



'A + Vl 



9m, 



2 1/ 1 + 2/2 
462 



mQ cos -j/^ 



2-2/^-2/2' 



^2 = 3mo(|/i - 7/2) 



2> C0S7 + 2a cos?/'2 



2a6 



2 , 2 2\'^'^°S'?/'2 — COS 7 



3"^o(yl -2/2 
"^0(2/? - y'2 



2ab 



2,2 2\'^si'^^2 + sin 7 



5(4 = mg(2 - yi - ^2) sinV'i , 

o 2/ 2 2Nsin7 - 2asinV2+ 
^6 = 3mo(yi - 1/2) 



(67) 



and 7 = '01 — '02 = tan~^ hmit of yi,|/2 — > 1 the parameters in Eq. ([6 

behave as p — > 3mQ,g — )■ 6mQ/a^,r — )• and (yfj — 0. So, as expected, the matrix in 

Eq. ( 165!) gives the TBM values. Similarly Eq. ( l25l) . we have three mixing angles (^^ 5 ^2; ^3)) 
three phases {(p'^, (f)2, (ps), and the three mass-squared eigenvalues. In turn, this mixing 
matrix f/^ becomes one of the mixing matrix composing the PMNS matrix. To see how the 
neutrino mass matrix given by Eq. fl62l) can lead to deviations from their TBM values, we 
first introduce three small quantities e^, {i = 1, 2, 3), which are responsible for the deviations 



of the 9i from their TBM values: 



l + ei, 9^2 = ^2, = sin-i (^-L 

Then the mixing matrix Ui, up to order Ei can be written as 

/ V2-e3 l+e3V2 



£3 ■ 



(68) 



V3 

(l+£i+£3V2)e'"^3 



V6 



+ 



£26 



V3 



V3 

\f2+e\\f2-e-i ege 
V6 



£26 



v 



(l-£i-£3V2)e~'''^i+'''3) £2e-''»2 (V2-£3-V2ei)e~''^i 



£26'' 



v/3 



V6 



V2 
l+£l 
V2 



(69) 
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Now, the straightforward calculation with the general parametrization of 11^ in Eq. ( 125|) 
leads to the expressions for the masses and mixing parameters 

Im [C] sin 0^ - Re [C] cos 0^ 



tan6';' 



tan 



tan 29^ 



lm[B] cos((/)^ - 0^) + Re[B] sm{(j)'( - 0^ 



arg(G) 



A, -A 



arg {c^^C + e'^h^^B) 



A2 — Ai 



arg(Z) 



(70) 



where = cos6'J', = sin6'j^, and 



= Acf - |c^C + e*'^i<5|sin2^2^ + A3S 



u2 
2 ' 



A2 = Fcf - \G\sm2ei + Ks'( 



u2 



= Kef + \G\ sm2ei + Fs'i 



u2 



And the squared-mass eigenvalues are given by 



K-F 



|G|cos2^i^ 



(71) 



m. 



\ „u2 \ „u2 \ „v2 \ „u2 

A1C3 — A2S3 2 ^2C3 — AiSj 



COS 2^^ 



m, 



cos 2^^ 



A3 + \c';C + e"f'h';B\ tan ■ (72) 



As is well-known, because of the observed hierarchy Am\^^^^ ^ ^"^Li = "^2 ~ ''^1 > 0' ^^"^ 
the requirement of a Mikheyev-Smirnov-Wolfenstein resonance for solar neutrinos, there are 
two possible neutrino mass spectra: (i) the normal mass hierarchy (NMH) rrii < m2 < nis, 
and (ii) the inverted mass hierarchy (IMH) < mi < 7712- The solar and atmospheric 
mass-squared differences are given by 

,2 



m 



1^2 



m 



2\Z\ 
sin 2^^ 




m: , for NMH 



^2 



m- , for IMH 



(73) 



which are constrained by the neutrino oscillation experimental results. We will discuss it 
numerically in the next section. 

Plugging Eqs. f l37|l and f l69|) into Eq. (160|) . the PMNS matrix is recast to 

( f/,n-^AV^'^f/.21 f/.12-^AV<^^t/.22 f/.13-^£AV*3[/^23 \ 
UyZl Uu-i2 

where the phase 0| is given as Eq. 0341) . From Eq. fl74p , the neutrino mixing parameters can 
be displayed in terms of the standard parametrization 



PMNS 



/13 



V 



Pu , (74) 



8| as 



sin^^ 
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e2 
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1 - 




t/e3 


2 



sin^ 6. 
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1 - 




U 


e3 


2 



sin 9 



13 



\U, 



e3| 



(75) 
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FIG. 1: The reactor mixing angle ^13 versus the ratio of third-to-second generation neutrino Yukawa 
couplings 1/1/1/2 (left plot) and the parameter k = y'^v^/M (right plot). The (red) crosses and 
(blue) dots represent the results for the inverted and the normal mass hierarchy, respectively. The 
horizontal dotted (solid) lines in both plots indicate the upper and lower bounds on ^13 given in 
Eq. ([H) at the 3c7 level (best-fit value). 

Leptonic CP violation at low energies can be detected through the neutrino oscillations 
which are sensitive to the Dirac CP-phase, but insensitive to the Majorana CP-phases in 



PMNS 



22|: the Jarlskog invariant [20| is defined as 



Jcp = -lm\Ul^Ue-i,UTiU*^] = ^ sin 26^12 sin 26^13 sin 26^23 cos 6^13 sin 6cp , (76) 

o 

where Uaj is an element of the PMNS matrix in Eq. fl74p . with a = e, corresponding to 
the lepton flavors and j = 1, 2, 3 corresponding to the light neutrino mass eigenstates. And 
by manipulation of Eqs. fl6Tl) and fl76l) one can easily obtain the Dirac CP phase : 

.CP = -arg I """""" - . (77) 



S12S23 



As expected, since the contributions of ¥[ in Eq. fl37|l to the PMNS matrix are negligible, 
i.e., its effects C(A^), we will consider f/pMNS — in numerical analysis. 



IV. NUMERICAL ANALYSIS 



Now we perform a numerical analysis for neutrinos using the linear algebra tools in 



Ref. 



23|. 
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FIG. 2: The reactor mixing angle ^13 versus the phase 93. The (red) crosses and (blue) dots 
represent the results for the inverted and normal mass hierarchy, respectively. The vertical dotted 
(solid) lines indicate the upper and lower bounds on ^13 given in Eq. ([1]) at the 3a level (best-fit 
value) . 

The mass matrices niD and Mr in Eq. fl62|) contains seven parameters: 
yi^,Vis,, M,yi,y2, The first three (7/3, M and f$) lead to the overall neutrino scale 

parameter mg. The next four {yi,y2, K,,ip) give rise to the deviations from TBM as well 
as the CP phases and corrections to the masse eigenvalues [see Eq. (IMl) ]. Since we have 
a relation v^/A ~ in the charged fermion sector, for the cutoff scale A = 10^^ GeV we 
take M = 10^^ GeV and vq, = 172\^ GeV, for simplicity, as inputs. Since the neutrino 
masses are sensitive to the combination mo = w||?/3^|/(6M), other choices of M and f$ give 
identical results. Then the parameters mQ,yi,y2, n,(p can be determined from the exper- 
imental results of three mixing angles, ^12,^13,^23; and the two mass squared differences, 
Amg^i, Am\^j^. In addition, the CP phases Sep, y^i,2 can be predicted after determining the 
model parameters. (Here, we will not discuss the Majorana CP phases v^i,2-) 

Using the formulae for the neutrino mixing angles and masses and our values of M,Vip, 
we obtain the following allowed regions of the unknown model parameters: for the normal 
mass hierarchy (NMH), 

0.21 <K< 0.86 , 1.04 < 2/1 < 1.37 , 1.04 < y^ < 1.39 , 

96° < V9 < 121° and 239° < < 257° , 1.1 < mo x lO^^ieV] < 3.6 ; (78) 
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FIG. 3: The behaviors of ^23 and 9i2 in terms of ^13. The red crosses and the blue dots represent 
results for the inverted mass hierarchy and the normal mass hierarchy, respectively. The dotted 
vertical lines represent the experimental bounds of Eq. ([1]) at Sir's. 

for the inverted mass hierarchy (IMH), 

0.81 <K<2, 0.91 < < 1.09 , 0.91 < 1/2 < 1.08 , 

116° < ^ < 248° , 1.6 < mo X lO^^ieV] < 2.2 . (79) 

Note that here we have used the Scr experimental bounds on ^12, ^23, Amg^j, Am^^j^ in 
Eq. ([1]), except for ^13 < 12° for which we use the values in Eqs. fl78|79p . For these parameter 
regions, we investigate how mixing parameters do behave for the NMH and IMH. In Figs. [1]- 
El the data points represented by blue dots and red crosses indicate results for the NMH 
and IMH, respectively. The left- hand- side plot in Fig. [1] shows how the mixing angle 6'i3 
depends on the ratio yi/y2 = Hi / 112 third- and second-generation neutrino Yukawa 

couplings; the right- hand- side plot shows how 6*13 depends on the parameter k = y'^v^/M. 
Fig. [2] shows the mixing angle ^13 as a function of the phase (p of y'^v^/M. As can be seen 
in Figs. [UIJl only normal mass hierarchy is permitted within 3a experimental bounds. And 
we see that the measured value of 6*13 from the Daya Bay and RENO experiments can be 
achieved at 3a's for 0.92 < y,/y2 < 1, 1< yi/y2 < 1.06, 0.2 < k < 0.7, 95° < ^ < 115° and 
245° < ^ < 255° for NMH. For IMH, in Figs. [T]l2] the value of ^13 reaches at most 3°, which 
is excluded by the measurements of ^13. 

Fig. [3] shows how the values of ^13 depend on the mixing angles ^23 and 612. As can be 
seen in the left plot of Fig. |3l the behavior of 6'23 in terms of the measured values of 6is at 
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FIG. 4: The Jarlskog invariant Jcp versus the reactor angle ^13 (left plot), and the Dirac CP phase 
5cp versus 623 (right plot). The (red) crosses and (blue) dots represent the results for the inverted 
and normal mass hierarchy, respectively. The vertical dotted (solid) lines indicate the upper and 
lower bounds on ^13 given in Eq. ([T]) at the 3a level (best-fit value). 



3cr's for the NMH is different than for the IMH. As aheady mentioned, the IMH is excluded 
by the measured values of ^13 in Fig. [31 For the NMH we see that the measured values of 
6'i3 can be achieved for 49.5° < 6*23 < 54.8° and 35.8° < 6*23 < 40.5°, with large deviations 
from maximality, which are favored at la by the experimental bounds as can be seen in 
Eq. ([T]). Future precise measurements of 6'23, whether 6*23 — )■ 45° or 16*23 ~ 45° | — )■ 5°, will 
either exclude or favor our model. From the right plot of Fig. [31 we see that the predictions 
for ^13 do not strongly depend on 6*12 in the allowed region. 

To see how the parameters are correlated with low-energy CP violation observables mea- 
surable through neutrino oscillations, we consider the leptonic CP violation parameter de- 
fined by the Jarlskog invariant in Eq. ( [76l) which can be expressed in terms of the elements 



of the matrix h = M^Ml 



22|: 



J ^ Im{/li2/j.23^3l} 

'"^ Am2iAm3]^Am|2 



(80) 



The behavior of Jcp as a function of ^13 is plotted on the left plot of Fig. [H We see that the 
value of Jcp lies in the ranges 0.006 ~ 0.03 and —0.03 ~ —0.004 (NMH) for the measured 
value of ^13 at 3cr's. Also, in our model we have 

27m^ 



Im{/li2/l23^3l} 



(81) 
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FIG. 5: Plots of \mee\ as a function of ^13 and mughtcst- The red crosses and the blue dots represent 
results for the inverted and the normal mass hierarchy, respectively. The vertical dashed lines show 
the experimental bounds of Eq. ([1]) at Sir's. 

in which { } stands for a comphcated lengthy function of yi, 1/2, a, b, ipi and 1^2- Clearly, 

Eq. (181 p indicates that in the limit of ?/2 — ^ Vi or sin'?/;2 the leptonic CP violation Jqp 
goes to zero. When y2 7^ yi, i.e. for the IMH case, Jcp could go to zero as smip2 of Eq. flSTj) 
[see, Eq. fl52l) and Fig. [2]. In the case of the NMH, Jcp has nonzero values for the measured 
range of 613 while Jcp goes to zero for ^13 — t- 0, which corresponds to y2 ^ yi- The right 
plot of Fig. m shows the behavior of the Dirac CP phase 6cp [see Eq. fl77j) ] as a function of 
^23, where the values of 6cp lie in the ranges 0° < 6cp < 60°, 130° < 6cp < 180°, 180° < 
S^p < 240° and 310° < 6cp < 360° for the NMH (for the IMH, 6cp can vary over a wide 
range, but which is excluded by the measured values of 613). Interestingly, for the best-fit 
values of 623 the values of 6cp are predicted as the one around 10°, 170°, 190°, 350° for NMH. 
So, future precise measurements of 6*23 will provide more information on 6c p- 

Moreover, we can straightforwardly obtain the effective neutrino mass |mee| that charac- 
terizes the amplitude for neutrinoless double beta decay [25|: 



irrif. 



PMNSjei"^* 



(82) 



where f/pMNS is given in a good approximation as Eq. fIBI?]) . The left and right plots in 
Fig. O show the behavior of the effective neutrino mass |mee| in terms of 613 and the lightest 
neutrino mass, respectively. In the left plot of Fig. [5l for the measured values of ^13 at 3(t's, 
the effective neutrino mass \mee\ can be in the range 0.04 < |mee|[eV] < 0.11 for NMH. The 
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right plot of Fig. |5] shows | function of mughtcst) where Slightest = fTT-i for the NMH 

and miightest = ^3 for the IMH. Our model predicts that the effective mass |mee| is within 
the sensitivity about 10~^ eV of planned neutrinoless double-beta decay experiments [26 1. 



V. CONCLUSION 



Under SU{2)l x U{1)y gauge symmetry, we have proposed a new model of leptons and 
quarks based on the discrete flavor symmetry T', the double covering of /I4. Here we impose 
that all Yukawa couplings be of order one, which implies that the hierarchies of charged 
fermion masses and the mildness of neutrino masses are responsible for six types of Higgs 
scalars. In addition to the gauge and flavor symmetries, in order to simplify our model and 
to remove the unwanted Yukawa terms appearing in the Lagrangian we have introduced 
a continuous global symmetry U{l)x which can not be gauged. After spontaneous U{l)x 
breaking, to avoid Goldstone bosons it has to be explicitly broken down to a subgroup 
Z2. After spontaneous breaking of flavor symmetry, with the constraint of renormalizability 
in the Lagrangian, the leptons have mg = and the quarks have CKM mixing angles 

= 13°, = 0° and 6*^3 = 0°. Thus, certain effective dimension-5 operators driven by the 
gauge-singlet and T'-triplet x fisld are introduced as an equal footing, which induce me 7^ 
and lead the quark mixing matrix to the CKM one in form. On the other hand, the neutrino 
Lagrangian still keeps renormalizability. We have assumed that there is a cutoff scale A, 
above which there exists unknown physics. 

We have shown numerical analysis in the lepton sector of our model, where only normal 
mass hierarchy is permitted within 3a experimental bounds with the prediction of both large 
deviations from maximality in the atmospheric mixing angle 623 and the measured values of 
reactor angle. So, future precise measurements of 623, whether 6^23 — )■ 45° or |023 — 45°| -> 5°, 
would either exclude or favor our model. Together with it, our model has made predictions 
both for the Dirac CP phase 0° < 5cp < 60°, 130° < 5cp < 180°, 180° < 5cp < 240° and 
310° < 6cp < 360°, which is almost compatible with the global analysis in la experimental 
bounds. Moreover, we have shown the effective mass \mee\ measurable in neutrinoless double 
beta decay to be in the range 0.04 < |mee| [eV^] < 0.11 for the normal hierarchy, which could 
be tested in near future neutrino experiments. 
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Appendix A: The Higgs potential 





X 




7] 



H 

G 



y = 



y = L 



V = L/2 

FIG. 6: Fifth dimension and locations of scalar and fermion fields. 

In this Appendix, as an example, we present our Higgs potential and its minimization, 
as well as our prescription for effecting the stability of the vacuum alignment. We solve the 



vacuum alignment problem by extending the model into a spatial extra dimension y [24 1. 
We assume that each field lives on a 4D brane either at y = or at ?/ = L/2 or at ?/ = L, 
as shown in Fig. |6l The heavy neutrino masses and the neutrino Yukawa interactions arise 
from local operators at the branes y = and y = L/2, while the charged fermion masses 
are realized by non-local effects involving both branes y = 0, y = L/2 and y = L. We 
impose that all the parameters appearing in the Lagrangian are assumed to be real. Once 
the scalars $, \E', x, H, G acquire complex VEVs at the different branes, the CP symmetry 
can be spontaneously broken. A rigorous explanation of this possibility is beyond the scope 
of this paper. 

The most general scalar potential {d < 5) for the Higgs fields $, H, \E', G, rj and x invariant 
under SU{2)l x U{1)y x T' and obeying the conditions in the previous paragraph, is then 
given by 

V = Vy=o + K,=f + Vy=L, (Al) 

where 



Vy=0 



y=L 



V{^) + V{^) + V{m) 
V{x) + V{r,) + V{xri) 
V{H) + V{G) + V{HG) 



(A2) 
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and 



V{^) = /i|(<l>t$)i + Af(<l>t<l>)i(<l>t$)i + Af($^<l>)i.($^$)i// 

+ Af($t$)3.($t$)3^ + Af($t$)3^($t$)3^ + Af{($t$)3^($t$)3^ + /,.c.} , (A3) 

V{^) = + Af + X^{¥^)r{^^^)i" + A*(^t^)3.(^%)3. 

+ AnM/tM/)3„(M/tM;)3^ + A^{(M/tM;)3^(M/tM;)3^ + h.c.} , (A4) 
^($^) = Af*(<l>t$)i(^t^)^ + Af*($t$),,($t^)^„ + Af*($t$)3^(^t^)3^ 
+ \ri^^^hS^^^)s^ + \r{i<^^<^)sX¥^)s^ + h.c.} 
+ A**{($t$)3„(*t^)3. + h.c.} + A^ {($^vl;)i($tvl>)i + h.c.} 
+ Xr + h.c.} + Af* {($^^)3.($^^)3. + h.c.} 

+ Afo* {(<^>^^)3„(<f^^)3„ + h.c.} + Aff {($tM/)3^($tM;)3^ + h.c.} 

+ Af2*($tv[/)i(v[/t$)^ + Af3* {($tvI/)i,(M/t$),„ + h.c.} + Aff (<l>tvl>)3^(x]>t$)3^ 

+ Af5'($^*)3.(^^$)3. + Af6*{($^*)3.(^^$)3. + /^.C.} , (A5) 

V{H) = ^^^UH^H)^ + X^iH^HUH^H)^ + X^ {{H^ HUH^ H)s + /i.e.} , (A6) 

= z/,2^(G'tG')i + Xf{G^G)^{G^G)l + Af {(G'tG')3(GtG)3 + /i.e.} , (AT) 

V{HG) = Xf^{H^H)i{G^G)i + X^^iH^G)i{G^H)^ + Xf^{{G^H)i{G^H)i + h.c.} 

+ Xf'^iiH^HUG^Gh + h.c.} + Af ^{(i/tG)3(GtiJ)3 + h.c.} 

+ Af ^{(Gt/J)3(Gt/7)3 + h.c.} , (A8) 

^(X) = /iJ{(xx)i + (xV)i} + rnJ(xX*)i + A^{to)i(xx)i + (x*X*)i(x*X*)i} 
+ A^ + ix*x)i"{x*X*)i'} 

+ Af {(xx)3.(xx)3. + ix*x*hAx*x*)3s} + Af (X*X)3. {(XX)3. + ix*X*)3s} 

+ A^ {(x*x)3„(xx)3. + {xx*hAx*x*hs} 

+ {x(xx)3. + x*{x*xl3s} + {x(x*x*)3. + x*{xx)3s} , 

+ f + {x*X*)i{x*X*xli} + - + %...} , (A9) 



In V{H) and y(G) the terms ifij^iW H)i and i/i^(GtG)i are expanded as i^\{HlH2 - hIHi) and 
iliQ(G\G2 — GjGi), respectively. 
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V{7i) = f^liv^v) + y^iv^vf , (AlO) 

v{xv) = ^r{v^v){{xx)i + {x*x*)i} 

+ j^iv^v){{xxx)i + {x*x*x*)i} + j^{v^v){ix*xx)i + {xx*x*)i} ■ (All) 

Here, fi^, fin, A^G; Z^??) A^x' "^x' ^i^ ^i^ Cl..i3 ^1^2 have a mass dimension, whereas Af 5, 
A* 5, X^2y ^1,2^ '^i,...,4? -^2,35 '^f,*.,i6' -^i!.^,4 "^^^ ^h dimensionless. And in V{x) 
" • • • " denotes dimension-5 operators composed of all possible combinations of x fields. 



1. Minimization of the neutral scalar potential 

After the breaking of the flavor and electroweak symmetry, in the neutral Higgs sector, 
in order to find minimum configuration of the Higgs potential, we in general let 

m = I , V I ' (^.) = I 1 V J , = f / . 

{H,) = [ ^ V I , (Gfc) =1 V I , (x,) = v^^e^^ , (A12) 

with j = 1 — 3, = 1, 2, where vq,^,v^,.,Vn,VH,,,VG,,Vy^^ are real and positive, and 7^, Cj; 
PkyCTkjfj are physically meaningful phases. Note that we can set ^ = without loss of 
generality because 9 does not have physical meanings [see, Eqs. ( IA9p - (lAlip ]. First, at the 
brane y = the vacuum configuration for $ and \& is obtained by vanishing of the derivative 
of V with respect to each component of the scalar fields $j and "^j. Then, we have six 
minimization equations for VEVs and six equations for phases. From those equations, we 
can get ^ 

- - 18aU4»1^„, (*.) = (*3) = 0, (A13) 



2{9Xf + AXI J 

where W = 9Af * + 4Af * + 9Af2* + 4Af * + 2(9Af" * + 4Af *) cos 2(7 - C) with 7 = 71, C = Ci, 
and v<j, and v^, are real. With the vacuum alignments of $, ^' fields, Eq. f lA13p . minimal 
condition with respect to 7j, Q are given as 



dV 



dV 



2„,2 



* * (9A** + 4AD sin 2(7 - C) = , (A14) 



9 



Of course, there are trivial solutions — 0,u*^ — 0. We have neglected them. 
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where 7 = 71, C = Ci? and = = is automatically satisfied. 

Second, at the brane y = L/2 the vacuum configuration for x ^iid rj is obtained by 
vanishing of the derivative of V with respect to each component of the scalar fields Xj a-^id 
77. For simplicity, we consider only the renormalizable terms in V{x) ^i-iid V{xv)- Then, we 
have seven minimization equations for four VEVs and three phases. From those equations, 
we can get ^ 



V 



+ 2(/i2 +i;2a>^'')cos2(/p 



Xj 



-1 + ^ 



(A15) 



12{(A^ + A^) cos4(^ + A^ cos 2^} 

where v-^ is real, and y^i = (/)2 = V^s = is used. With the vacuum alignment of x fields, 
Eq. flAlSp . minimal condition with respect to (/:>j is given for i{>\ = i{)2 = </?3 as 

1 dV 



|t;2A>^'' + fil + 3vl (a^ + 4(A^ + A^) cos 2ip,^ } sin 2ip, = , (A16) 



with i = 1,2, 3. And, requiring vanishing of the derivative of V with respect to 77, 



1 dV 



2dri^ 



<X>='"x 



•^A'' + ^ + 3A*>^^;2^os2<^ 



M^'v ■ 2 



<r;">=j;^ 

and, we obtain the VEV of t] for (x) = t>^e*'^(l, 1, 1), 

2 _ -fil-QvlX^''cos2ip 



(A17) 



2A^ 



(A18) 



Finally, at the brane y = L the vacuum configuration for H and G is obtained by vanishing 
of the derivative of V with respect to each component of the scalar fields Hi and Gj. Then, 
we have eight minimization equations for four VEVs and four phases. From those equations, 
we can get 



1 dV <G°>=VG^ 



2dHf 



1 dV 



2 dm 



<GO>=VG, 



<Hp=VH^ 



t;^,,|2i;|^Af (cos2pi2 - 1) - Af^i;^^ 



- fGif G2 (2Af *^ sin (712 sin pi2 - Af^ cos(cri2 + P12)) } = ) 
= fiy2|2t;^^Af (cos2pi2 - 1) - AfS^^j + i;^! |ph sinpi2 

- vg.vg^ (2Af^sinai2sinpi2 - A^^ cos((7i2 + P12)) } = , (A19) 



^ There are trivial solutions v^. — Q. We have neglected them. 
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1 dV 



2 



1 dV 



2 



<H'>>=v„^ 



t'Gi|2fG2-^f (cos2ai2 - 1) - AfSI^I + t^Gaj/ic sin 0-12 
vh^vh^ (2Af'^ sin (T12 sin P12 - Af'^cos(cri2 + P12)) } = ' 
^Gajs^^Gi^f (cos2cri2 - 1) - Af | + t;^ j/ic sm CT12 
WHi^^i?2 (2Af ^ sin cri2 sin pi2 - A^^ cos((Ti2 + P12)) } = • (^20) 



where pi2 = Pi — P2 and (7i2 = (Ti — 0-2. And, we obtain the VEVs of H and G: 

sinpi2 - vg^vg^Y 



VG2 



2t;^,,Af (1 - cos2pi2) + X^^'vIJvh, 

2vH,\U^ - cos2pi2) + X^^vlJvH, 
/i^ sin 0-12 - VHiVHjy 

2vG,\f{l - COS 2^12) + \2^vjjJvG, 

p^sin(Ti2 - vh^vh^Y 



(A21) 



2t;G,Af (cos2(Ti2 - 1) - X^^'vIJvg, ' 
where y = 2Af^'^ sin cri2 sin pi2 — A^*-^ cos((Ti2 + Pi2)- With the vacuum ahgnment of H and 



G fields, Eqs. ( ]A19P and (lA20p . minimal condition with respect to pjjCTj is given as 



1 dV 


1 dV 






~ 2dp2 


= vh.vh^^^ - 






+ vg^vg2 (2A 


1 dV 


1 dV 




2dai 


~ 2da2 





+ vh,vh2 (2Af^sinpi2Cosai2 + Af^sin(pi2 + ^12)) } = • (A23) 



Appendix B 



In Eq. (IT5]) the components m^- are given by 



i-i 



-VH2 + , 



'12 



21 



m'22 = Yl 







— i 




2 






— i 


I 2 



(Bl) 
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In Eq. (fT6|) the components mfj are given by 

mil = -YiVH, - Y^^-^vh, + 17 {^vh, - ^-^^H^ 

mi, = -Y.^VH, + 17 (^^'^^2 - ^^H,) - Y^ (^^vh, + ^^^^2) 



mi. 



"^32 



nV. - + (^^H, - ^^H.) , mi, = (f - f ) • (B2) 

And, in Eq. f lT7|) the components m*^- are given by 

mil = Y: (^^-j^vg, - - Y: (^^-J^^'g^ + ^vg,^ 

= Y:^-^vg, + Y: (^IvG, - i^^c.) , m\, = v,^f 
mil = Y:'-vg^ + Y: (-^^Gi + ^^^2) 

< = -Y:\vg, + Y: (^t'G. - IvG,^ , m*3 = (f + f ) , 

"^31 = -y:^^^vg, + y: {jvg^ — 1 , 

ml, = Y: Qf;o. - ^^c.) - (i^^c, + 1^,,,) , m*3 = (f - f ) (B3) 
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